home *** CD-ROM | disk | FTP | other *** search
/ Language/OS - Multiplatform Resource Library / LANGUAGE OS.iso / dsp / dr.bub / 96000.lha / 96000 / appb / b125.asm < prev    next >
Assembly Source File  |  1992-04-28  |  4KB  |  57 lines
  1. ; This program was originally published in the Motorola DSP96002 Users Manual
  2. ; and is provided under a DISCLAIMER OF WARRANTY available from Motorola DSP
  3. ; Operation, 6501 William Cannon Drive West, Austin, Texas 78735-8598.  For
  4. ; more information, refer to the DSP96002 Users Manual, Appendix B, DSP
  5. ; Benchmarks.
  6. ;
  7. ; B.1.25    Argument Reduction  
  8. ;Argument reduction (AR) is the problem of having a desired floating  point number range and an ar-
  9. ;gument that is outside of the range.  The  argument is placed inside of the desired range by adding or 
  10. ;subtracting  multiples of the desired number range.  Of course, adding and  subtracting multiples of a 
  11. ;number is inherently slow and requires  infinite precision. Some simple methods can be used with 
  12. ;some  assumptions on the precision of the data and relative argument sizes.  
  13. ;The following program performs AR when the desired range is arbitrary  and the input value is arbi-
  14. ;trary. This may be used to reduce an angle to  the range of -pi to pi.  
  15. ;The following variables are defined: 
  16. ; 
  17. ;           rmin = range minimum value, -pi in this example 
  18. ;           rmax = range maximum value, pi in this example 
  19. ;           range = rmax-rmin, 2*pi in this example 
  20. ;           o_range = 1.0/range 
  21. ;
  22. ;Assume the input is in d0.  
  23. ; 
  24. ;rmin      equ  -3.14159 
  25. range     equ  2*3.14159 
  26. o_range   equ  1.0/range 
  27.  ;                                                          Program    ICycles 
  28.  ;                                                           Words 
  29.  
  30.  move             #range,d7.s   ;load desired range 
  31.  move             #rmin,d2.s    ;load range min 
  32.  move             #o_range,d3.s ;load reciprocal of range 
  33.  
  34.  fadd.s  d2,d0                  ;adjust to rmin               1     1 
  35.  fmpy.s  d0,d3,d0               ;scale the input              1     1 
  36.  floor   d0,d1                  ;get integer part             1     1 
  37.  fsub.s  d1,d0                  ;get fractional part          1     1 
  38.  fmpy.s  d7,d0d0                ;spread out fraction to range 1     1 
  39.  fadd.s  d2,d0                  ;adjust to rmin               1     1 
  40.   ;                                                           ---   --- 
  41.   ;                                                   Totals:  6     6 
  42.  
  43. ;The output is in d0. Note that the constant initialization is not  included in the benchmark because it 
  44. ;does not need to be executed every  time argument reduction is desired and is therefore application  
  45. ;dependent.  
  46. ;If the desired range begins at zero (i.e. the desired range is zero to  two pi), then the first and last 
  47. ;fadd instructions can be deleted for a  four cycle argument reduction.  
  48. ;This is one possible method for AR and it is efficient.  This method  will not work when the argument 
  49. ;divided by the result range has no  fractional part (in the current precision).  This is obvious since it  
  50. ;is the fractional part that contains the information relating to how  far the scaled argument is in the 
  51. ;reduced range. The integer part tells  how many times the range has wrapped around.  Typically, a 
  52. ;good  programmer will keep the argument to a few multiples of the desired  range. In most practical 
  53. ;applications, the argument may exceed the  desired range by several integral values.  In this case, the 
  54. ;presented  algorithms work very well.  After the final reduced argument has been  obtained, any incre-
  55. ;ments should be made from the reduced argument to  prevent eventual overflow and maintain maxi-
  56. ;mum precision.  
  57.